Doping and temperature dependence of superfiuid weights for high T c cuprates 
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Using the improved U(l) and SU(2) slave-boson approaches of the t-J Hamiltonian [Phys. Rev. B 
64, 052501 (2001)] that we developed recently, we report the doping and temperature dependence of 
the superfiuid weight. It is shown that at low hole doping concentration x and at low temperatures T 
there exists a propensity of a linear decrease of the superfiuid weight n s /m* with temperature, and 
a tendency of doping independence in the linearly decreasing slope of -^-(x,T) vs. T in agreement 
with the experimentally observed relation ^(x,T) = -^(x,0) — aT with a, a constant. It is also 
demonstrated that both T c and n s /m* increase with hole doping concentration x in the underdoped 
region, reaches a saturation (maximum) at a hole doping above optimal doping and decreases beyond 
the saturation point in the overdoped region. Such a reflex (decreasing) behavior of T c and n s /m* 
is attributed to the weakening of coupling between the spin(spinon pair order) and charge(holon 
pair order) degrees of freedom in the overdoped region. All of these findings are in agreement with 
jj.SK measurements. 
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Presently there exist two outstanding physical prob- 
lems of superfiuid weights to be resolved. One is the un- 
derstanding of the so-called 'boomerang' (reflex) behav- 
ior and the other, the doping independence of linear de- 
crease of superfiuid weight zrf^(x,T) vs. temperature T. 
More than a decade ago the transverse field muon-spin- 
relaxation(/i-SR) measurements of the magnetic penetra- 
tion depth A in high T c copper oxide superconductors re- 
vealed the 'boomerang' behavior of the superconducting 
temperature T c and the superfiuid weight n s /m*. that 
is, a linear increase of both T c and n s /m* in the under- 
doped region, a saturation of T c at optimal doping and a 
decrease (reflex) of both T c and n s /m* in the overdoped 
region as the hole doping concentration increases [1-3]. 
However, theoretical explanations of this boomerang (re- 
flex) behavior are still incomplete. Most recently, mag- 
netic penetration depth(A) measurements [4-6] showed a 
doping independence in the linearly decreasing slope of 
the superfiuid weight with increasing temperature obey- 
ing the relations ^fi-(x, T) = ^{x, 0) — aT in the under- 
doped region where x is the hole concentration and a, a 
constant. In these experiments [4-6], the normalized su- 
perfiuid weight shows a universal behavior as a function 



of scaled temperature, that is, 



n s / m* (x,T) 



1-0(T/T e ). 



n s /m* (a;,0) 

This universal behavior has been interpreted by the phase 
fluctuation model [7] , or the quasiparticlc model based on 
'the Landau- Fermi liquid theory [8-10]. However, apart 
.from these phenomenological models the theoretical ex- 
planation of this universal behavior from a microscopic 
theory is largely unaccomplished [11]. 

■ 

There have been numerous efforts to theoretically ex- 
plain the microscopic origin for this doping independent 
propensity of the decreasing slope of the observed super- 
fluid weight with increasing temperature [12-14]. How- 
ever, until now the problem still remains unexplained 
[11]. Earlier U(l) slave-boson theories [12,13,15-17] of 
the t — J Hamiltonian predicted that at small hole dop- 
ing the superfiuid weight ^{x,T) as a function of tem- 



perature is strongly doping (x 2 ) dependent. This is in 
direct conflict with the /iSK measurements [4-6]. Lee 
and Wen [12,18] proposed that the SU(2) theory which 
incorporates the low energy phase fluctuations of order 
parameters may resolve this problem. Recently D.-H. 
Lee [13] presented a U(l) slave-boson theory concerned 
with the low energy fluctuations of order parameters and 
the excitations of massless Dirac fermions at the d-wave 
nodal point. This theory, also, showed the x 2 depen- 
dence of at finite temperature. Most recently Wang et 
al. [14] also showed that the recently proposed d-density 
wave theory [19] fails to predict the doping independence 
of the linearly decreasing slope of superfiuid density vs. 
temperature. 

Earlier, we presented an improved U(l) and SU(2) 
slave-boson approach of the t-J Hamiltonian which differs 
from other slave-boson approaches to the t-J Hamiltonian 
in that coupling between the charge and spin degrees of 
freedom is manifested in the slave-boson representation 
of the Heisenberg term of JJ2 <i j > {Si ■ Sj — \niTij), by 
naturally admitting the charge degree of freedom exhib- 
ited in the intersite charge coupling energy term of jUiTij 
[20]. In this theory [20] holons (charges) are subject to 
attractive interaction as a result of coupling to the spinon 
singlet pairs. As a result, holons form holon pairs and the 
holon-pair bosons undergo bose condensation [21] at low 
temperatures. This is the major difference of the present 
holon-pair boson theory as compared to other earlier the- 
ories which pay attention to singe-holon bose condensa- 
tion [12,13,15-18]. Our theory was shown to predict the 
arch shaped superconducting (holon-pair bose condensa- 
tion) temperature T c as a function of doped hole con- 
centration x [20], in agreement with the experimentally 
observed phase diagram. It is, thus, of great interest to 
employ this theory for checking its predictability of the 
boomerang and both the doping and temperature depen- 
dence of the superfiuid weight. 

Since the predicted results of various physical proper- 
ties are qualitatively indistinguishable between the U(l) 
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and SU(2) slave-boson theories, below we present only 
the SU(2) theory for generality. In the SU(2) slave- 
boson representation [18], the electron operator is given 

( /i 

by c Q = -^h^ipa with a = 1,2, where ipi = 



U'2 



h 
~fl 



and h 



are respectively the dou- 



blets of spinon and holon annihilation operators in the 
SU(2) theory. The SU(2) slave-boson representation of 
the t-J Hamiltonian shows [20] 



H = 



"2 E [tfLuM 3 bu - blb 2] ) + h.c. 

<i,3> 

Hhihj - fufvWljbx + blb 2j ) + h.c. 
-f E {^ + h\h i ){l + h]h j )x 

<i,j> 

i 

E [^(flfl + h \M + ^f\hihi + blfiu) 



(1) 



where A) 



(1),(2),(3) 



are the real Lagrangian multipliers to 



enforce the local single occupancy constraint in the SU(2) 
slave-boson representation [18]. It is noted that the in- 
clusion of the charge-charge interaction —jTiirij in the 
Hcisenbcrg Hamiltonian results in the coupling between 
the spin (spinon) and charge (holon) degrees of freedom. 
It is easy to show that under inverse SU(2) transforma- 
tion of the above expression we recover the original form 
of the U(l) slave-boson representation given in the liter- 
ature [20]. 

After relevant Hubbard-Stratonovich transformations 
[20] , the holon action is obtained to be 

S b (A, X , Af, A B ,h) = f dr f V ht(Ti, r)(d T - fi)h(r t , r 
Jo 1 l 

4 E i4i 2 (Ei a ^i 2 +* 2 ) 



<hi> 



a, 



2 ^ 

<i,j> 



4 EO 

<i,j> 



[e^h\r t ,T)U^h{v 3 , T ) + c.c. 



|Af ? | 2 / 1 T(r l ,r)A-(/ l T (r3 , r) )- +c . c . 



(2) 



n ;m (w,q) =U b l ^' A) (uj,q) 

E E 

a 1 ,a 2 =e,a,!3l\ m '=x,y 



and the spinon action, 



r/3 r 

5 / (x,A/,V) = J o dr[Y,^{^r)d T ^{v i ,T) 



+ 



J(l 



E (i4i 2 + §^i 



<i,j> 



-{l-xf (^ t (ri,r)E^(r J -,r) + C . C .) 



(3) 



<i,j> 



Here various symbol definitions are as follows. /i(r i; r) = 
h(r if T) \ . g the Sl] ^ doublct of holon ficld; and 

/i(J"»,t) 




, the SU(2) doublet of spinon 



A is the electromagnetic (EM) field. \]} ■ = 



-A/.* 



A/. 

Xij 



and U/j = 



-2A 



f* 



~~ Xij 



1,3 



are re- 



spectively the order parameter matrices of hopping(xij) 
and spinon pairing(A^). Here Xi,i+l — r\e x<Xl ^ Ti ^ cos(#° + 
9i(ri)) and A f ii+l = rje 113 '^ sin(0° + 6» ; (r,)), with 77 = 
v/|x| 2 + |A/| 2 and 0° = itan^ 1 ^(thc sign +(-) is for 
the ij link parallel to x (y)). 07 (r^) represents the phase 
fluctuations of the hopping order parameter, /3;(rj), the 
phase fluctuations of the spinon pairing order parame- 
ter, and 6i(ri), the relative phase fluctuations between 
the amplitudes of hopping and spinon pairing order pa- 



rameters. A B 4 



A b A b 
A b A b 



is the matrix of the 



holon pairing order parameter. 

After integration over the holon and spinon fields, we 
obtain the total free energy, 



F(A) 



1 



In / D x DA f DA l 



-0(F\A, x ,A f A B )+F f (X-A 1 )) 



(4) 



where F b (A, x , A', A B ) = -I In/ £>/ ie -s b (A,x,A',A B ,h) 
is the holon free energy and Ff(x,Af) 
— j] In J Dipe~ s, ( x ' Af the spinon free energy. 

To compute the EM current response function, we first 
obtain the total free energy by integrating out the three 
phase fields a, (3 and 9 above. We obtain a formula for 
the EM current response function, up to the second or- 
der, 



q) [nV, q) + n V q)] v n S™ A) 



(5) 
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where III 6 (a;, q) + 11^ (u, q)l 2 / is the inverse ma- 



trix element of 



n £' a V,q) + nf/V a V,q) 



Here 



a ^ anc ^ ^-fm 1 Q ^ represent the isospin current re- 



sponse functions of holons and spinons respectively to 
gauge fields a 1 , a 2 = 9, a, or /3. Similarly, IL^' A \(j, q) 
is the EM current response function of holon to the EM 



field, n 



b(a\A) 
Im 



is the isospin current response function 



of holon to both gauge fields a and A. It is noted that 
the SU(2) isospin current of holon(spinon) is defined as 



Mf) 

Ja.l 

3i = -/3 



-P g^o, if> , and analogously for the 'EM' current, 
In the SU(2) slave-boson theory, the re- 



8F 

8 A, ■ 



sponse function n^/ 4 '" ^ of the holon isospin current to 
the EM field vanishes owing to the contribution of the &2- 
boson in the static and long-wavelength limit [18] . There- 
fore the superfluid weight of the total system is given only 
by the holon current response function, 



Us 

m* 



= 4limn^»(c = 



(w = 0,q). 



(6) 



Here n^'^ (oj, q) is computed from the use of the usual 
linear response theory for the holon action Eq.(2) [22]. 

In Fig. 1(a), with the use of our slave-boson theory [20] 
we show the temperature dependence of superfluid weight 
for a wide range of doping concentrations covering both 
underdoping and overdoping. The slope of each curve 
represents the variation of superfluid weight with tem- 
perature which is not in precise agreement with observa- 
tions [4-6], overestimating the superfluid weight ^{x, 0) 
at T = OK and further ^fi-(x,T) reaches zero rather 
abruptly compared to the observation which showed a 
slower drop [5]. We would like to point out that a close 
look at Fig. 1 of the above experimental paper [5] re- 
veals a linear slope behavior only in the lower temper- 
ature region and a quick non-linear drop of superfluid 
weight at the higher temperature region. Such a change 
of slope from linearity is qualitatively in agreement with 
our results although our theory predicted a faster drop in 
^r(T). The continuous drop in the predicted superfluid 
weight indicates a second order phase transition. We be- 
lieve that improvement can be achieved if we can intro- 
duce the t term in our t-J Hamiltonian which is known 
to cause a good agreement with the observed dispersion 
energy by showing realistic cold and hot spots. The con- 
tribution of the hot spot in the Brillouin zone becomes 
increasingly important with increasing temperature up 
to T c to allow classical phase fluctuations. 

As temperature increases, the predicted superfluid 
weight shows a tendency of linear decrease in T with 
nearly identical slopes particularly at low temperatures 
in the underdoped region. To rigorously check such a 
propensity of doping independence in the slope of ^ vs. 
T, we computed differences in the superfluid weight at 
T = and T ^ 0, that is, %r(x,T) - %fi:{x,0) at each 
x. The results arc shown in Fig. 1(b). A decrease 



of the superfluid weight with temperature at all dop- 
ings is predicted although the temperature dependence 
of the superfluid weight does not show a clear linear- 
ity at all temperatures as shown in Fig. 1 (b). There 
exists a tendency of linear decrease with a nearly iden- 
tical slope particularly in the underdoped region only at 
low temperature. Although the temperature dependence 
of 2*r(x,T) does not show a global linearity by yield- 
ing a good fit to the empirical relation of ^(x,T) = 
^(x,0) - a(x)T - a(x)T 2 - 0(T 3 ), the linear slope 
a(x) is found to be doping in dependent only in the un- 
derdoped region at low temperature as is clearly shown 
in Fig. 1 (b). As shown in the figure it is found that 
the variation of the slope, that is, the superfluid weight 
n s /m* vs. T/t at low temperatures does not appreciably 
change in the underdoped region (x = 0.04, 0.07 and 0.1). 
On the other hand, in the overdoped region (x = 0.16, 
0.18 and 0.2) the slope of the superfluid weight with tem- 
perature is no longer doping-independent quite unlike the 
underdoped case, again, in agreement with observation. 
The predicted optimal doping is x Q = 0.13. It is noted 
that there exists a crossing behavior between the un- 
derdoped and overdoped lines of superfluid weights vs. 
temperature as shown in Fig. 1(a). This prediction is 
consistent with measurements [4-6]. 

In Fig. 1(c) we display the the normalized superfluid 
weight ( " a /m*(z'o) ) as a f unc ti° n °f th e scaled tempera- 
ture (T/T c ) for various hole concentrations. As shown 
in the figure, the normalized superfluid weight begins 
to decrease linearly at low temperature and drops more 
rapidly at higher temperature, by showing a convex line 
shape in qualitative agreement with experiments [4-6]. 
The decrease of superfluid weight with temperature oc- 
curs as a consequence of diminishing spin singlet pairing 
order with increasing temperature. This is readily under- 
stood from the holon pairing term (the last term in Eq. 
(2)) which reveals the interplay of the spin (spinon) sin- 
glet pairing order with the holon pairing order; the spin 
singlet pairing order is predicted to decrease with temper- 
ature. Thus the decrease of the superfluid weight with in- 
creasing temperature occurs by breaking the holon pairs 
(hole pairs) in connection with the diminishing singlet 
pairing order at higher temperatures. A linear decrease 
of superfluid weight is predicted at low temperature and 
a rapid decrease at higher temperature, displaying a con- 
vex line shape. This trend of the convex line shape for all 
doping cases is in complete agreement with observation 
[4-6]. Although the predicted /? value in the relation of 

ra S /m*(x'o) = 1 ~~ fi{T/T c ) shows linear decrease particu- 
larly in the lower temperature region in agreement with 
observation, its (3 value scattered around 0.055 fails to 
quantitatively agree with the empirical value of (3 « 0.5 
independent of doping in the underdoped region [6] . Such 
discrepancy indicates that the inclusion of the next near- 
est hopping (t ) term is important. This is because the 
Fermi surface segment grows from the nodal point ini- 
tially at half-filling and closes near (it, 0) as hole doping 



3 



reaches optimal doping and beyond as is well-known from 
the angle resolved photoemission spectroscopy measure- 
ments [24]. In the future the inclusion of the t term is 
necessary to take a good account of classical phase fluctu- 
ations in the region of the hot spot at high temperatures 
near and at T c 

According to the present theory the superconductivity 
is accomplished by the condensation of the Cooper pairs 
of d-wave symmetry as composites of the holon pairs of 
s-wave symmetry and the spinon pairs of d-wave symme- 
try. Thus the Cooper pairs can be regarded as composite 
objects resulting from such coupling of the spin (spinon 
pairing) and charge (holon-pairing) degrees of freedom. 

To clearly show the effect of coupling between the 
spin(spinon) and chargc(holon) degrees of freedom, we 
calculated the superfluid weight by varying the antifer- 
romagnetic spin-spin coupling, that is, the Heiscnbcrg 
coupling J. In Fig. 2, we show the J dependence of su- 
perfluid weight at x = 0.1 for three different choices of 
J/t = 0.2, 0.3 and 0.4. As J increases, a decreasing ten- 
dency in the slope of the superfluid weight vs. tempera- 
ture is predicted. For stronger Heisenberg coupling J and 
thus the spinon pairing strength (order), it is more dif- 
ficult to thermally break the holon-pairs due to a larger 
holon-spinon coupling. Such coupling effect is manifested 
in the last term of Eq.(2). This results in the more slowly 
decreasing slope of the superfluid weight with increasing 
temperature for large J, as shown in the figure. 

In Fig. 3, we show the predicted boomerang behavior 
in the plane of T c vs. ^(x, T — > 0), by showing a linear 
relationship between the two in the underdoped region 
and the 'reflex' behavior in the overdoped region based 
on the U(l) theory. By the reflex behavior we mean the 
decrease of both T c and as the hole doping concen- 
tration increases in the overdoped region. This predicted 
trend is consistent with the measurements [1-3] of the 
muon-spin-relaxation rates a (a oc n s /m* cx T c ). The 
pseudogap temperature is predicted to rapidly decrease 
in the overdoped region in agreement with observation. 
The predicted pseudogap caused by the appearance of the 
spinon pairing order decreases in the overdoped re- 
gion and this, in turn, causes a decrease in holon pairing 
A fc and the coupling effect between the charge (holon- 
pair) and spin (spinon pair) degrees of freedom in the 
overdoped region thus a decrease in both T c and 
This coupling between A* and A b is exhibited in the last 
term of Eq.(2). Such coupling appears in exactly identi- 
cal manner in both the U(l) and SU(2) theory. Thus the 
boomerang behavior of the SU(2) theory is qualitatively 
the same as that of the U(l) theory. 

In summary, we have found the following qualita- 
tively salient features for the doping and temperature 
dependences of the superfluid weight. We demonstrated 
the nearly doping independence in the slope of super- 
fluid weight vs. temperature in the underdoped re- 
gion, thus satisfying the experimentally observed relation 
^ix.T) = £k(x,0) - aT with a, a constant. On the 
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other hand, in the overdoped region the doping depen- 
dence in the slope of the superfluid weight with temper- 
ature is correctly predicted in agreement with observa- 
tion. From the present study we found a tendency of the 
boomerang behavior in the locus of T c with the varia- 
tion of hole doping and superfluid weight in the plane of 
^r(x,T — > 0) vs. T c . This is in qualitative agreement 
with the ^i-SR experiments. The decreasing (boomerang) 
behavior of T c and ^ in the overdoped region is at- 
tributed to the weakened coupling of holon pairs to the 
spin pairing order. 
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FIG. 1. (a) Temperature dependence of the superfluid 
weight for underdoped(a; = 0.04,0.07,0.1), optimal 

doping(a; = 0.13) and overdoped(:r = 0.16,0.18,0.2) rates with 
J/t — 0.2. (b) The difference of superfluid weights between 
T / OA" and T = OA", that is, ^_(T) - ^(T = 0). Each 
line represents a fitting to the computed superfluid weight, (c) 
The normalized superfluid weight „'^»|ioj as a function of 
scaled temperature T/T c for same doping concentrations and J 
value as (a). 
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FIG. 2. The difference of superfluid weights between T / OK 
and T = OK, (T) - ^ (T = 0) for a underdoped hole concen- 
tration x = 0.1 with J/t = 0.2, 0.3 and 0.4. Each line represents 
a fitting to the computed superfluid weight. 
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FIG. 3. The superfluid weight ^(x, T) at T = 0.001* (equiv- 
alent to T ~ 5-R" with the use of * = 0.44eV[23]) vs. supercon- 
ducting temperature T c with the choice of J/t = 0.2 based on 
the U(l) theory. The open box represents hole concentration 
starting from x = 0.01 to x = 0.22 and the arrow denotes the 
direction of increasing doping rate. The predicted optimal dop- 
ing rate is x D = 0.07 with the U(l) theory. The SU(2) theory 
predicts qualitatively the same boomerang behavior. 
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